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Abstract 

We study the hadroproduction of two neutral Higgs bosons in the minimal su- 
persymmetric extension of the standard model (MSSM), which provides a handle on 
the trilinear Higgs couplings. We include the contributions from quark-antiquark 
annihilation at the tree level and those from gluon-gluon fusion, which proceeds 
via quark and squark loops. We list compact results for the tree- level partonic 
cross sections and the squark loop amplitudes, and we confirm previous results for 
the quark loop amplitudes. We quantitatively analyze the hadronic cross sections 
at the CERN Large Hadron Collider assuming a favorable supergravity-inspired 
MSSM scenario. 

PACS numbers: 12.60.Fr, 12.60.Jv, 13.85.-t 



1 Introduction 



The search for Higgs bosons will be among the prime tasks of the CERN Large Hadron Col- 
lider (LHC) [|I|. While the standard model (SM) of elementary-particle physics contains 
one complex Higgs doublet, from which one neutral CP-even Higgs boson H emerges 
in the physical particle spectrum after the spontaneous breakdown of the electroweak 
symmetry, the Higgs sector of the minimal supersymmetric extension of the SM (MSSM) 
consists of a two-Higgs-doublet model (2HDM) and accommodates five physical Higgs 
bosons: the neutral CP-even and bosons, the neutral CP-odd boson, and the 
charged iP'^-boson pair. At the tree level, the MSSM Higgs sector has two free parameters, 
which are usually taken to be the mass m^o of the boson and the ratio tan (3 = V2/V1 
of the vacuum expectation values of the two Higgs doublets. 

In the SM, the trilinear and quartic Higgs-boson couplings are proportional to the 
square of the Higgs-boson mass, and they are uniquely fixed once the latter is known. In 
the MSSM, the various Higgs-boson self-couplings are determined by the gauge couplings 
multiplied with trigonometric factors that depend on a, the mixing angle that rotates 
the weak CP-even Higgs eigenstates into the mass eigenstates and H^, and /? [see 
Eq. (Q]. 

At the LHC, the trilinear Higgs-boson couplings may be probed by studying the inclu- 
sive hadroproduction of Higgs-boson pairs. These couplings enter the stage via the Feyn- 
man diagrams where a virtual neutral Higgs boson is produced in the s channel through 
qq annihilation or gg fusion and in turn decays into a pair of neutral or charged Higgs 
bosons. Specifically, the partonic subprocesses include qq — > HH and gg —>■ HH in the 
SM, and qq 0i02, H+R- and gg 0i02, H+R- in the MSSM, where 0i = H^, A^. 
At the tree level, there are two mechanisms of qq annihilation. On the one hand, it can 
proceed via a Z boson (Drell-Yan process) if an appropriate Higgs-Higgs-Z coupling ex- 
ists. On the other hand, the Higgs bosons can be radiated off the g-quark line if the 
relevant Yukawa couplings are sufficiently strong. In the SM, qq annihilation is greatly 
suppressed due to the absence of a HHZ coupling and the smallness of the Hqq cou- 
plings for the active quarks contained inside the proton, q = u,d, s, c, b. In the MSSM, 
however, we have h^A^Z and H^A^Z couplings at the tree level [see Eq. ( |A.2| )], and the 
(pibb couplings are generally strong if tan/3 is large [see Eq. ( |A.3|) ]. In the SM, gg fusion is 
mediated via heavy-quark loops. In the MSSM, there are additional contributions from 
squark loops. The SM case was studied in Ref. [Q. As for H^H~ pair production, qq 
annihilation was investigated in Refs. the quark loop contribution to gg fusion in 

Refs. , and the squark loop one in Refs. . As for the pair production of neutral 

Higgs bosons in the MSSM, qq annihilation via a Z boson was analyzed in Ref. |§| and 
the quark and squark loop contributions to gg fusion in Refs. [Q and [|T0| , respectively. In 



Ref. [g[, also the QCD corrections to qq annihilation via a Z boson and to gg fusion via 
an infinitely heavy top quark were considered. The processes of associated production of 
a neutral Higgs-boson pair with a dijet (in addition to the remnant jets), an intermediate 
boson, or another neutral Higgs boson were found to have cross sections that are greatly 
suppressed, by more than an order of magnitude, compared to those of the respective 



2 



processes without those additional final-state particles []Tl] . 

The purpose of this paper is to reanalyze the pair production of neutral Higgs bosons 
in the MSSM, both via qq annihilation and gg fusion. In the case of qq annihilation, 
we also allow for q = b. This makes it necessary to include a new class of diagrams 
involving 6-quark Yukawa couplings, which are depicted in Fig. ^(a) and the second and 
third lines of Fig. |l|(b). These come in addition to the Drell-Yan diagram, shown in the 
first line of Fig. |l](b), which is already present for q = u,d, s, c. For final states with an 
even number of CP-even Higgs bosons, 

(j)^ = h'^h\ h^H°, H^H^, A^A^, bb annihilation 
is the only production mechanism at the tree level. As we shall see, its cross section is 
comparable to — and in certain areas of the MSSM parameter space even in excess of 
— the one of gg fusion. For final states with an odd number of CP-even Higgs bosons, 
0102 = hPA'^,H^A'^, bb annihilation significantly enhances the Drell-Yan contribution for 
large values of tan f3. As for gg fusion, we reproduce the analytical and numerical results 
for the quark loop contributions of Ref. [0. Apart from obvious typographical errors, we 
also agree with the formulas for the squark loop contributions listed in Ref. [|10|. However, 
we find their numerical size to be considerably smaller than what was found in Ref. pi| . 

As for bb annihilation, it should be noted that the treatment of bottom as an active 
fiavor inside the colliding hadrons leads to an effective description, which comprises con- 
tributions from the higher-order subprocesses gb — >■ 0i026, gb 0i02&, and gg 010266. 
If all these subprocesses are to be explicitly included along with bb 0i02, then it is 
necessary to employ a judiciously subtracted parton density function (PDF) for the b 
quark in order to avoid double counting [|12]. The evaluation of bb — » 0i02 with an un- 
subtracted 6-quark PDF is expected to slightly overestimate the true cross section . 
For simplicity, we shall nevertheless adopt this effective approach in our analysis, keeping 
in mind that a QCD-correction factor below unity is to be applied. 

In order to reduce the number of unknown supersymmetric input parameters, we 
adopt a scenario where the MSSM is embedded in a grand unified theory (GUT) involving 
supergravity (SUGRA) |T^. The MSSM thus constrained is characterized by the following 
parameters at the GUT scale, which come in addition to tan (3 and m^o: the universal 
scalar mass mo, the universal gaugino mass mi/2, the trilinear Higgs-sfermion coupling A, 
the bilinear Higgs coupling B, and the Higgs-higgsino mass parameter /i. Notice that m^o 
is then not an independent parameter anymore, but it is fixed through the renormalization 
group equation. The number of parameters can be further reduced by making additional 
assumptions. Unification of the r-lepton and 6-quark Yukawa couplings at the GUT scale 
leads to a correlation between rrit and tan f3. Furthermore, if the electroweak symmetry 
is broken radiatively, then B and fi are determined up to the sign of /i. Finally, it turns 
out that the MSSM parameters are nearly independent of the value of A, as long as 
|A| ^500 GeV at the GUT scale. 

This paper is organized as follows. In Sec. ^ we list analytic results for the tree-level 
cross sections of qq —>■ 0i02, including the Yukawa-enhanced contributions for q = b, and 
the squark loop contributions to the gg — * 0i02 amplitudes in the MSSM. The relevant 
MSSM coupling constants and the squark loop form factors are relegated to Appendices A 
and B, respectively. In Sec. |^, we present quantitative predictions for the inclusive cross 
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section of pp 0i02 + X at the LHC adopting a favorable SUGRA-inspired MSSM 
scenario. Sec. | contains our conclusions. 



2 Analytic Results 

In this section, we present the tree-level cross sections of the partonic subprocesses qq — > 
0102; where 0j = H'^, A^, and the transition (T) matrix elements of gg — >■ 0i02 arising 
from squark triangle and box diagrams. 

We work in the parton model of QCD with nf = 5 active quark flavors, which we take 
to be massless. However, we retain the 6-quark Yukawa couplings at their finite values, 
in order not to suppress possibly sizeable contributions. We adopt the MSSM Feynman 
rules from Ref. 0. In Appendix A, we list the trilinear self-couplings of the neutral Higgs 
bosons and their couplings to gauge bosons and quarks. For each quark flavor q there 
is a corresponding squark flavor q, which comes in two mass eigenstates i = 1,2. The 
masses rriq^ of the squarks and their trilinear couplings to the and bosons are listed 
in Eq. (A. 5) of Ref. and Eq. (A. 2) of Ref. p[, respectively. Their trilinear couplings 



to the boson and their quartic couplings to the H^, and A^ bosons may be found 
in Appendix A. 

Considering the generic partonic subprocess ab 4'i4>2, we denote the four-momenta 
of the incoming partons, a and b, and the outgoing Higgs bosons, (pi and 02, by pa, Pb, 
Pi, and p2, respectively, and define the partonic Mandelstam variables as s = {pa + PbY, 
t = {pa — Pi)^, and u = {pb — Pi)^ ■ The on-shell conditions read = = and pf = hi, 
where hi denotes the square of the 0i-boson mass. Four-momentum conservation implies 
that s + 1 + u = hi + h2. Furthermore, we have sp|. = tu — hih2, where pt is the absolute 
value of transverse momentum common to 0i and 02 in the center-of-mass (cm.) frame. 

The tree-level diagrams for bb 0i02, with 0i02 = h^h^,h^H^,H^H^,A^A^ and 
0102 = h^A^^H^A^, are depicted in Figs. |I](a) and (b), respectively. The cross sections 
for the first class of partonic subprocesses may be generically written as 

where the prefactor accounts for identical-particle symmetrization if 0i = 02, is 
Fermi's constant, my/ is the VT-boson mass. 



•S* = g<t,i<t,2hoghObb'Phois) + g<j,:^<j,2HogHObb'PHois), 



T± — g<j)^bbg<t,2bb ( 7 =t - ) ■ (2) 



Here, 



Vx{s) = ^— — (3) 



iln Ref. 0, TUq^ is called M^^. 
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is the propagator function of particle X, with mass mx and total decay width Tx- For 
the second class of partonic subprocesses, we have 

^{bb M2) = [\P\'+ pU\V\' + \A~ T+n] , (4) 

where 

P = 9(l,i4,2AogAObbPAo{s), 

V = 2g^^^.^zVzbb'Pz{s), 

A = 2g^^^.^zazbb'Pz{s). (5) 

Here, vztb = -(/&- 2slQb) / {2cui) and azbb = -/f,/(2c„), with cl = 1 - si = my^/m\, 
are the vector and axial- vector couplings of the b quark, with weak isospin /j, = —1/2 
and electric charge Qb = —1/3, to the Z boson. As for h'^A^ and H^A^ production, there 
are also sizeable contributions from qq annihilation via a Z boson for the quarks of the 
first and second generations, q = u,d, s, c. The corresponding Drell-Yan cross sections 
are obtained from Eq. (|^) by putting P = T+ = and substituting b ^ q. The resulting 
expression agrees with Eq. (36) of Ref. [Q. The full tree-level cross sections are then 
obtained by complementing the 66-initiated cross sections of Eq. with the Drell-Yan 
cross sections for q = u,d, s, c. 

The one-loop diagrams for gg —>■ 0i02, with 0i02 = h^h^,h^H^,H^H^,A^A^ and 
0102 = h^A^,H^A^, are depicted in Figs. |^(a) and (b), respectively. As for the quark 
loops, our analytical results fully agree with those listed in Ref. 0, and there is no 
need to repeat them here. For the partonic subprocesses of class two, the squark loop 
contributions are zero 0]. This may be understood as follows, (i) The ggg^g^ and gzgtgj 
couplings are linear in the squark four-momenta, while the gggg^gj couplings are momentum 
independent. Thus, the diagrams in the third line of Fig. |^(b) each vanish upon adding 
their counterparts with the loop- momentum flows reversed, (ii) The ggg^g^ , Qggqigj , QhOqigj , 
and gHOgicjj couplings are symmetric in i and j, while the gA'^Mj coupling is antisymmetric. 
Thus, the diagrams in the last line of Fig. ^(b) vanish upon summation over i and j. For 
the partonic subprocesses of class one, the T-matrix elements corresponding to the squark 
triangle and box diagrams are found to be 

= ^el{p.)el{pb)ATF^, 

V2 TT 

fn = ^elip.)e:ipb) {a^Fu + AfGo) , (6) 

respectively, where asifJ^r) is the strong-coupling constant at renormalization scale fir, 
^liiPa) is the polarization four- vector of gluon a and similarly for gluon b, it is summed 
over the color index c = 1, . . . , 8, 

2 

^1 —9 PbPai 

s 

=9'' + \ i-p'bP: + ^^p'ip: + ^pCk + p^k) , (7) 
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and the form factors Fa, Fn, and Gn are listed in Appendix B. Due to Bose symmetry, 
7a and To are invariant under the simultaneous replacements fi ^ u and Pa ^ Pb- 
Consequently, Fa, Fn, and Gn are symmetric in t and u. Our analytic results for the 



squark loop contributions agree with those given in Eqs. (8)-(10) of Ref. [|T0|, which are 
expressed in terms of helicity amplitudes.^ 

The parton- level cross section of gg — * 0i02 including both quark and squark contri- 
butions is then given by 



da . 



1 Gj^ajifir 
1 + 6^,^, 256(27r): 



+ 



Gn — 



2m 



w 



G, 



Q=t,b 



{fa + h 



where the generalized couplings and and the form factors F^, Fn, and Gn may 
be found in Eq. (16)-(18) and Appendix A of Ref. 0, respectively. 

The kinematics of the inclusive reaction AB GD + X, where A and B are hadrons, 
which are taken to be massless, and G and D are massive particles, is described in Sec. II of 
Ref. |T^. Its double-differential cross section d'^a/dy dp^, where y and pt are the rapidity 



and transverse momentum of particle G in the cm. frame of the hadronic collision, may 
be evaluated from Eq. (2.1) of Ref. HI 



3 Phenomenological Implications 

We are now in a position to explore the phenomenological implications of our results. 
The SM input parameters for our numerical analysis are taken to be Gp = 1.16639 x 
10-5 GeV-2, mw = 80.419 GeV, mz = 91.1882 GeV, rrit = 174.3 GeV , and rub = 4.6 GeV 
]1§. We adopt the lowest-order (LO) proton PDF set CTEQ5L fT^. We evaluate as{fir 



from the LO formula |jT6| with Uf = 5 quark flavors and asymptotic scale parameter 
^QCD — 146 MeV []T^. We identify the renormalization and factorization scales with the 
0102 invariant mass ^/s, M = /i^ = ^/s. We vary tan f3 and m^o in the ranges 3 < tan f3 < 
38 ~ nit /nib and 90 GeV < niAo < 1 TeV, respectively. As for the GUT parameters, we 
choose mi/2 = 150 GeV, A = 0, and /i < 0, and tune mo so as to be consistent with 
the desired value of itiao- All other MSSM parameters are then determined according to 
the SUGRA-inspired scenario as implemented in the program package SUSPECT |T^ . 



For the typical example of tan (3 = ?> and m^o = 300 GeV, the residual masses and total 
decay widths of the 0i bosons are m^iO = 90 GeV, nijjo = 306 GeV, P^o = 3 MeV, 
Tho = 186 MeV, and P^o = 72 MeV, and the squark masses are m^^ = nic^ = 412 GeV, 
niu^ = = 422 GeV, mj^ = ms^ = 413 GeV, mj^ = ms^ = 428 GeV, = 317 GeV, 
mj^ = 443 GeV, m^^ = 384 GeV, and m^^ = 413 GeV, We do not impose the unification of 

^ There are two obvious typographical errors on the right-hand side of Eq. (lOe) in Ref. There 
should be an overall minus sign, and V/f^. -yq^q^ should be replaced by Vh^^. ^^q^qi ■ 
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the r-lepton and 6-quark Yukawa couphngs at the GUT scale, which would just constrain 
the allowed tan/? range without any visible effect on the results for these values of tan/3. 
We exclude solutions which do not comply with the present experimental lower mass 
bounds of the sfermions, charginos, neutralinos, and Higgs bosons |]19|. In our analysis, 
an s-channel resonance only occurs in the process pp — *• h^h^ + X if nifjo > 2mho. 

We now study the fully integrated cross sections of pp — > </)i</)2 + X at the LHC, with 
cm. energy = 14 TeV. Figures |-| refer to the cases (/)i</)2 = h^h°, h^H^, H^H^, 
h^A^, H^A^, respectively. In part (a) of each figure, the 771^0 dependence is studied for 
tan/3 = 3 and 30 while, in part (b), the tan/3 dependence is studied for m^o = 300 GeV. 
We note that the SUGRA-inspired MSSM with our choice of input parameters does not 
permit tan /? and m^o to be simultaneously small, due to the experimental lower bound 
on the selectron mass |T^. This explains why the curves for tan/3 = 3 in Figs. ^|-^(a) 



only start at rriy^o 240 GeV, while those for tan /? = 30 already start at m^o 90 GeV. 
On the other hand, tan/3 and m^o cannot be simultaneously large either, due to the 
experimental lower bounds on the chargino and neutralino masses [0. For this reason, 
the curves for tan/3 = 30 in Figs. |3-B|(a) already end at m^o ^ 560 GeV. Finally, the 



experimental m/jo lower bound |]I9[ enforces tan /3 ^ 3 if m^o = 300 GeV, which is reflected 
in Figs. |-|(b). 

In Figs. 0-^, the 66-annihilation contributions (dashed lines), which originate from 
Yukawa-enhanced amplitudes, and the total gg-fusion contributions (solid lines), corre- 
sponding to the coherent superposition of quark and squark loop amplitudes, are presented 
separately. For a comparison with future experimental data, they should be added. For 
comparison, also the gg-fusion contributions due to quark loops only (dotted lines) are 
shown. We first assess the relative importance of the 66-annihilation and (^(^-fusion contri- 
butions. In the case of h^h^ production, bb annihilation is more important than gg fusion 
for intermediate values of m^o, around 300 GeV, except at the edges of the allowed tan/3 
range, while it is greatly suppressed for large values of m^o, independent of tan/3 (see 
Fig. In the case of h^H^ production, bb annihilation dominates for intermediate to 
large values of rriAo and large values of tan/3, while it yields an insignificant contribution 
for small values of tan /3, independent of m^o (see Fig. ^). As for H^H^ and A^A^ produc- 
tion, bb annihilation is suppressed compared to gg fusion. For tan/3^8, the suppression 
factor is modest, ranging between 2 and 3, but it dramatically increases as tan/? becomes 
smaller (see Figs. |^ and In order to avoid confusion, we should mention that the 
66-annihilation contribution for tan/? = 3 is too small to be visible in Figs. ^-^(a). We 
now investigate the size of the supersymmetric corrections to gg fusion, i.e., the effect of 
including the squark loops. We observe that these corrections can be of either sign and 
reach a magnitude of up to 90%. Specifically, for h'^h^, h^H^, H^H^, and A^A^ produc- 
tion, they vary within the ranges -10% to +3%, -16% to -^32%, -32% to +24%, and 0% 
to +90%, respectively, for the values of uiao and tan /3 considered in Figs. 0-|^(a). The fact 
that they are relatively modest in most cases is characteristic for our SUGRA-inspired 
MSSM scenario. This is partly due to the destructive interference of quark and squark 
loop amplitudes and to the suppression of the latter by heavy-squark propagators. By 
contrast, in Ref. []TU|, the supersymmetric corrections were found reach values in excess 
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of 100. Since the authors of Ref. did not specify all their input parameters, we could 
not reproduce their numerical results for the squark loop contributions. 

In Figs. 1^ and the total gg-annihilation contributions (dashed lines), corresponding 
to the coherent superposition of Drell-Yan and Yukawa-enhanced amplitudes, and the gg- 
fusion contributions (solid lines), which now only receive contributions from quark loops, 
are given separately. For comparison, also the pure Drell-Yan contributions (dotted lines) 
are shown. Again, we first compare the total gg-annihilation contributions with the gg- 
fusion ones. In the case of h^A^ production, qq annihilation dominates for large values 
of tan/?, independent of m^o, while, at the lower end of the allowed tan/3 range, it is 
suppressed by a factor of 40 and more, depending on the value of m^o. On the other 
hand, in the case of H^A^ production, the gg-annihilation contribution always overshoots 
the gg-iusion one by at least one order of magnitude. We then examine the effect of 
including the Yukawa-enhanced amplitudes in the evaluation of the gg-annihilation cross 
section. In the case of h^A^ production, there is a dramatic enhancement for large values 
of tan/5, which may reach several orders of magnitude for large values of m^o. In the case 
of H^A^ production, there is also an enhancement for large values of tan /5, but it is much 
more moderate, less than a factor of three. It is interesting to observe that the Drell-Yan 
cross section of H^A^ production is fairly independent of tan (3 unless m^o is close to its 
lower bound. This may be understood by observing that sin(a — /3), which governs the 
H^A^Z coupling gH°A°Zi defined in Eq. (|A.2|) , is then always very close to —1. This is 
also apparent from Fig. 2 of Ref. . 

At this point, we should estimate the theoretical uncertainties in our predictions. As 
a typical example, we consider the cross section of pp — > h'^h'^ + X for tan /3 = 3 and 
m^o = 300 GeV. In order to obtain a hint on the size of the as-yet unknown next-to- 
leading-order (NLO) corrections, we define the renormalization and factorization scales 
as M = /ir = ^^/s and vary the scale parameter in the range 1/2 < ^ <2. The resulting 
variation in cross section amounts to ±8% in the case of hh annihilation and to ±11% in 
the case of gg fusion. At NLO, one also needs to specify a renormalization scheme for 
the definition of the 6-quark mass, which enters our analysis through the 6-quark Yukawa 
coupling. Our LO analysis is appropriate for the on-mass-shell scheme, which uses pole 
masses as basic parameters. The modified minimal-subtraction (MS) scheme |^ provides 
a popular alternative. For example, a pole mass of m;, = 4.6 GeV corresponds 
to an MS mass of ml^\jjr) = 2.7 GeV for the typical choice of renormalization scale 
fir = = 300 GeV. Recalling that the leading behaviour of Eq. (|l|) in nih is quadratic, 
switching to the MS scheme would thus, at first sight, lead to a suppression of the cross 
section by a factor of approximately 1/3. However, we must keep in mind that this 
reduction should be largely cancelled, up to terms that are formally beyond NLO, by a 
respective shift in the NLO correction. Another source of uncertainty is related to the 
choice of PDF's. In fact, there exist significant differences in the extraction of the 6-quark 
PDF among different PDF sets, which are related to the threshold treatment of the g ^ bb 
splitting, the choice of the 6-quark mass, the dependence of the evolution on the latter, 
etc. In the case of a LO analysis with Uf = 5 massless quark favours, which is considered 
here, the use of CTEQ5L fl^ together with the scale choice M = fir = \/s should be 
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appropriate in the sense that these issues can largely be bypassed. If we employ the LO 
PDF set by Martin, Roberts, Stirling, and Thorne (MRST) H, with A'§.^ = 132 MeV, 
then the 66-annihilation and gfgf-fusion cross sections increase by 10% and 4%, respectively, 
relative to their default values. 



4 Conclusions 

We analytically calculated the cross sections of the partonic subprocesses qq — > 0i02 and 
gg 0102, where (pi = h'^,H^,A^, to LO in the MSSM. We included the Drell-Yan 
and Yukawa-enhanced contributions to qq annihilation (see Fig. P and the quark and 
squark loop contributions to gg fusion (see Fig. ^. We listed our formulas for the qq- 
annihilation cross sections and the squark loop amplitudes, for which we found rather 
compact expressions. As for the quark loop contributions, we found complete agreement 
with Ref. 0. 

We then quantitatively investigated the inclusive cross sections of pp — >• 0i02 + X 
at the LHC adopting a favorable SUGRA-inspired MSSM scenario, varying the input 
parameters m^o and tan/?. The results are presented in Figs. ^j-^. We found that the 
Yukawa-enhanced gg-annihilation contribution, which had previously been neglected, can 
play a leading role, especially for h^h^ production if m^o is of order 300 GeV and for 
hPH^, hPA^, and H^A^ production if tan/3 is large. The supersymmetric corrections to gg 
fusion, which are present for h^h^ , h^H^, H^H^, and A^A^ production, can be of either 
sign and reach a magnitude of up to 90%. Our numerical results for these corrections 
disagree with those presented in Ref. [|1^]. For each process pp —>■ (j)i(j)2 + X, the combined 



cross section, i.e., the sum of the full gg-annihilation and gg-fusion contributions, varies 
by several orders of magnitude as the values of m^o and tan f3 are changed within their 
allowed ranges, and its maximum value is typically between 10^ fb~^ and 10^ fb~^. If we 
assume the integrated luminosity per year to be at its design value of L = 100 fb~^ for 
each of the two LHC experiments, ATLAS and CMS, then this translates into a maximum 
of 20.000 to 200.000 events per year for each of these signal processes. 

A comprehensive discussion of the background processes competing with the pp 
(f)i(f)2 + X signals at the LHC lies beyond the scope of our study. However, we should 
briefly mention them and quote the relevant literature. Without specifying the decay 
channels of the (pi bosons, one expects the major backgrounds to arise from the pair 
production of neutral gauge bosons, the associate production of a neutral gauge boson 
and a neutral Higgs boson, and the continuum production of the respective (pi(p2 decay 
products. Published signal-to-background analyses have concentrated on the 0102 
bbbb signals and their irreducible continuum backgrounds, due to the partonic subprocesses 
gg,qq — > bbbb, which are dominantly of pure QCD origin. It has been demonstrated that, 
after optimizing the acceptance cuts, the LHC experiments might discover a signal, with 
experimental significance in excess of 5, if tan /? ^ 3 or tan (3 ^ 50. 
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A Relevant Higgs and squark couplings 

In this appendix, we hst the trihnear self-couphngs of the h^, and bosons as well 
as their couplings to the Z boson and the t and b quarks. Furthermore, we collect the 
couplings of these Higgs bosons to the squarks g^, with q = t,b and i = 1,2, which are 
not contained in Appendix A of Ref. [^]. For convenience, we introduce the short-hand 
notations Sa = sin a, Cq, = cos a, S/j = sin/3, cp = cos/5, S2/3 = sin(2/5), C2/3 = cos(2/5), 
s± = sin(a ± (3), and c± = cos(a ± (3). 

The trilinear self-couplings of the h^, H^, and A° bosons are given by 0] 

3mz mz ,^ . 

ghPh^hO = C2aS+, ghOhORO = --^ [^S2aS+ - C2aC+), 

mz , . ^mz 

ghOHORO = {^S2aC+ + C2aS+), gH^H^HO = ^ C2aC+, 

mz mz , . 

ghOAOAO = C2/3S+, gHOAOAO = C2/3C+. (A.lj 

Their couplings to the Z boson are given by P| 

c_ s_ 
ghOA^z = 7^ — , gn^A^z = 7, — ■ (A-2) 

Their couplings to the t and h quarks are given by 

ghott = —7; ' ghObb 



mfSa mbCo 

gnott = —- , gmtb 



nit cot (3 mfetan/? 

gAOtt = — , gAObb = — • (A. 3) 

The missing couplings of these Higgs bosons to the squarks are given by 
f gA^hh gA0hi2 \ ^ mt{^i + At cot/?) / l\ 
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\9A0b2bi 9A%2b2 . 

I 9h0h»ilil 9h0h0ili2 
\9h0h0i2ii 9h0h0i2i2. 

( 9h0h0bibi 9h0h%ib2 
\9h0h0b2bi 9h0h%2b2 , 

I 9hOHHih 9h°H°hi2 

\9h0H0t2h 9 hO HO 1212, 

( 9h0H0bibi 9hOH%ib2 
\9hPH%2bi 9h°H%2b2 . 

9H0H0iiii 9mH0iii2 
V 9H0H0t2ti 9 HO HO 1212 . 

9H0H0bibi 9H0H0bib2 
,9H0H0b2bi 9H0H0b2b2, 

( 9A0A0hh 9A0A0ht2 

\9AOAOi2h 9A0A0t2t2, 

( 9A0A0bibi 9A0A0bib2 
\9A0A0b2bi 9A0AOb2b2. 



mi,{fj, + Ah tan /?) / 


l\ 






-10) 




/ C2a(^It-sl,Qt) 




AW? 

= M 


W 


— — 2 2" 




\ 
\ 






/ C2a{h~sljQb) 




= m'^ 








I " 








mls2a 
















/ S2a{h-stQb) 1 
o„2 "T 


mls2a 



C2aS^Qt 

2cl 



2cl 



2m^ 



M' 



x\ T 



S2asl,Qt 

2c2, 



r\ T 





C2a (-ft-S^Qt) 



2c2 





+ 



T 



= M' 


24 


2m^.2 




/ C2a(4-s^Qf,) 


9 9 


= M^ 


2r2 


2^14/ '^z? 




^ C2/3(/t-S^Qt) 


cot^ /3 


= M' 


24 


2m^ 


= M^ 


/ C2f3(Ib-siQb) 


tan^ /3 
2m^ 



2r2 











„ 2 2 


2r2 






° \ 


'^2^^wQt 


m^cot^/J 1 


2r2 


2m^ / 










2r2 


2m^ / 



(A^^)'^.(A.4) 



Here, denotes the mixing matrix which rotates the left- and right-handed squark 
fields, qL and g/j, into the mass eigenstates g^. Its definition may be found in Eq. (A.l) 
of Ref . P . Relations similar to Eq. ( |A.4| ) are valid for the squarks of the first and second 
generations, which are also included in our analysis. However, in these cases, we neglect 
terms which are suppressed by the smallness of the corresponding light-quark masses. 



B Squark loop form factors 

In this appendix, we express the squark triangle and box form factors. Fa, Fu, and Gn, 
in terms of the standard scalar three- and four-point functions, which we abbreviate as 
Q/fc(c) = Cq (a, 6, c, m|, m\,ml^ and A'ifcf (e, /) = ^0 (a, 6, c, d, e, /, m| , m|, , m?^ , m| 
respectively. The definitions of the latter may be found in Eq. (5) of Ref. . 
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We have 

2 



q 1=1 



~ 2 

= - 9<Pmqj9<t>2qjqiF2 (s, /ll, /l2, , 



~ 2 

= 9ct>mqj9<t>2qjqiF3 (s, t, hi, /l2, ^| , ^|) • (B-l) 

Here, we have introduced the auxihary functions 

Fi{s,ml)=2 + 4mlC^i{s), 



F, (s, t, hi, h2, m|, m|) = -tiC^^'it) - hC^jfit) + 2smlD^}t'\s, + s f + 



Pt I 2 



Fa t, /.I, h„ ml, m|.) = (t + m|) C^{s) + ^m? Cg^^l^) - - tt^Cjf (t) 

+ (t^ - hih,)C^^l^-{s) + [st^ - 2iit2m| + 2sml (m| - m|.)] 
X D^^:^^'\s, t) - 2stmlD^l^'\s, t) + sm| (p| + m| - m|.) 
xDSf^°(t,«) + (t^«), (B.2) 

where ti — t — hi. 
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b h 



b f 



-H.A'^ 



-H.,A^ 



b f 



• H. 



H. 



(a) 



(b) 



Figure 1: Tree-level Feynman diagrams for qq — > 0i02, with (a) 0i02 
hPhP, h°H\ H'^H^ and (b) 0102 = in the MSSM. 
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(a) (b) 

Figure 2: One-loop Feynman diagrams for gg — > 0i02, with (a) 0i02 = 

h^h'^,h^H^,H^H°,A°A^ and (b) 0i02 = due to virtual quarks and squarks 
in the MSSM. 



16 





(b) 

Figure 3: Total cross sections a (in fb) of pp — > h^h'^ + X via bb annihilation (dashed lines) 
and gg fusion (solid lines) at the LHC (a) as functions of m/^o for tan/9 = 3 (starting at 
m^o = 240 GeV) and 30 (starting at m^o = 90 GeV); and (b) as functions of tan/5 for 
m^o = 300 GeV. For comparison, also the quark loop contribution to gg fusion (dotted 
lines) is shown. 
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Figure 4: Total cross sections a (in fb) of hS^H^ + X via 66 annihilation (dashed 

lines) and gg fusion (solid lines) at the LHC (a) as functions of myio for tan /5 = 3 (starting 
at m^o = 240 GeV) and 30 (starting at m^o = 90 GeV); and (b) as functions of tan/3 for 
rriAO = 300 GeV. For comparison, also the quark loop contribution to gg fusion (dotted 
lines) is shown. 
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Figure 5: Total cross sections a (in fb) of H^H^ + X via 66 annihilation (dashed 

lines) and gg fusion (solid lines) at the LHC (a) as functions of myio for tan /5 = 3 (starting 
at m^o = 240 GeV) and 30 (starting at m^o = 90 GeV); and (b) as functions of tan/3 for 
rriAO = 300 GeV. For comparison, also the quark loop contribution to gg fusion (dotted 
lines) is shown. 
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(b) 

Figure 6: Total cross sections a (in fb) of — > + X via 66 annihilation (dashed 

lines) and gg fusion (solid lines) at the LHC (a) as functions of myio for tan /5 = 3 (starting 
at m^o = 240 GeV) and 30 (starting at m^o = 90 GeV); and (b) as functions of tan/3 for 
rriAO = 300 GeV. For comparison, also the quark loop contribution to gg fusion (dotted 
lines) is shown. 
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Figure 7: Total cross sections cr (in fb) of pp h^A'^ + X via qq annihilation (dashed 
lines) and gg fusion (solid lines) at the LHC (a) as functions of my^o for tan /? = 3 (starting 
at m^o = 240 GeV) and 30 (starting at m^o = 90 GeV); and (b) as functions of tan/3 
for TTiAo = 300 GeV. For comparison, also the Drell-Yan contribution to qq annihilation 
(dotted lines) is shown. 
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Figure 8: Total cross sections a (in fb) of pp — > H'^A^ + X via qq annihilation (dashed 
lines) and gg fusion (solid lines) at the LHC (a) as functions of m^o for tan (3 = 3 (starting 
at m^o = 240 GeV) and 30 (starting at m^o = 90 GeV); and (b) as functions of tan/3 
for m^o = 300 GeV. For comparison, also the Drell-Yan contribution to qq annihilation 
(dotted lines) is shown. 
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